RNA degradome analysis reveals DNE1 endoribonuclease is required for the turnover of diverse mRNA substrates in Arabidopsis.

In plants, cytoplasmic mRNA decay is critical for posttranscriptionally controlling gene expression and for maintaining cellular RNA homeostasis. Arabidopsis DCP1-ASSOCIATED NYN ENDORIBONUCLEASE 1 (DNE1) is a cytoplasmic mRNA decay factor that interacts with proteins involved in mRNA decapping and nonsense-mediated mRNA decay (NMD). There is limited information on the functional role of DNE1 in RNA turnover, and the identities of its endogenous targets are unknown. In this study, we utilized RNA degradome approaches to globally investigate DNE1 substrates. Monophosphorylated 5' ends, produced by DNE1, should accumulate in mutants lacking the cytoplasmic exoribonuclease XRN4, but be absent from DNE1 and XRN4 double mutants. In seedlings, we identified over 200 such transcripts, most of which reflect cleavage within coding regions. While most DNE1 targets were NMD-insensitive, some were upstream ORF (uORF)-containing and NMD-sensitive transcripts, indicating that this endoribonuclease is required for turnover of a diverse set of mRNAs. Transgenic plants expressing DNE1 cDNA with an active-site mutation in the endoribonuclease domain abolished the in planta cleavage of transcripts, demonstrating that DNE1 endoribonuclease activity is required for cleavage. Our work provides key insights into the identity of DNE1 substrates and enhances our understanding of DNE1-mediated mRNA decay.

Integrative Analysis of Gene Expression and miRNAs Reveal Biological Pathways Associated with Bud Paradormancy and Endodormancy in Grapevine.

Transition of grapevine buds from paradormancy to endodormancy is coordinated by changes in gene expression, phytohormones, transcription factors, and other molecular regulators, but the mechanisms involved in transcriptional and post-transcriptional regulation of dormancy stages are not well delineated. To identify potential regulatory targets, an integrative analysis of differential gene expression profiles and their inverse relationships with miRNA abundance was performed in paradormant (long day (LD) 15 h) or endodormant (short day (SD), 13 h) Vitis riparia buds. There were 400 up- and 936 downregulated differentially expressed genes in SD relative to LD budsGene set and gene ontology enrichment analysis indicated that hormone signaling and cell cycling genes were downregulated in SD relative to LD buds. miRNA abundance and inverse expression analyses of miRNA target genes indicated increased abundance of miRNAs that negatively regulate genes involved with cell cycle and meristem development in endodormant buds and miRNAs targeting starch metabolism related genes in paradormant buds. Analysis of interactions between abundant miRNAs and transcription factors identified a network with coinciding regulation of cell cycle and epigenetic regulation related genes in SD buds. This network provides evidence for cross regulation occurring between miRNA and transcription factors both upstream and downstream of MYB3R1.

Murine Oviductosomes (OVS) microRNA profiling during the estrous cycle: Delivery of OVS-borne microRNAs to sperm where miR-34c-5p localizes at the centrosome.

Oviductosomes (OVS) are nano-sized extracellular vesicles secreted in the oviductal luminal fluid by oviductal epithelial cells and known to be involved in sperm capacitation and fertility. Although they have been shown to transfer encapsulated proteins to sperm, cargo constituents other than proteins have not been identified. Using next-generation sequencing, we demonstrate that OVS are carriers of microRNAs (miRNAs), with 272 detected throughout the estrous cycle. Of the 50 most abundant, 6 (12%) and 2 (4%) were expressed at significantly higher levels (P < 0.05) at metestrus/diestrus and proestrus/estrus. RT-qPCR showed that selected miRNAs are present in oviductal epithelial cells in significantly (P < 0.05) lower abundance than in OVS, indicating selective miRNA packaging. The majority (64%) of the top 25 OVS miRNAs are present in sperm. These miRNAs' potential target list is enriched with transcription factors, transcription regulators, and protein kinases and there are several embryonic developmentally-related genes. Importantly, OVS can deliver to sperm miRNAs, including miR-34c-5p which is essential for the first cleavage and is solely sperm-derived in the zygote. Z-stack of confocal images of sperm co-incubated with OVS loaded with labeled miRNAs showed the intracellular location of the delivered miRNAs. Interestingly, individual miRNAs were predominantly localized in specific head compartments, with miR-34c-5p being highly concentrated at the centrosome where it is known to function. These results, for the first time, demonstrate OVS' ability to contribute to the sperm's miRNA repertoire (an important role for solely sperm-derived zygotic miRNAs) and the physiological relevance of an OVS-borne miRNA that is delivered to sperm.

Analysis of Brachypodium miRNA targets: evidence for diverse control during stress and conservation in bioenergy crops.

BACKGROUND: Since the proposal of Brachypodium distachyon as a model for the grasses, over 500 Bdi-miRNAs have been annotated in miRBase making Brachypodium second in number only to rice. Other monocots, such as switchgrass, are completely absent from the miRBase database. While a significant number of miRNAs have been identified which are highly conserved across plants, little research has been done with respect to the conservation of miRNA targets. Plant responses to abiotic stresses are regulated by diverse pathways many of which involve miRNAs; however, it can be difficult to identify miRNA guided gene regulation when the miRNA is not the primary regulator of the target mRNA. RESULTS: To investigate miRNA target conservation and stress response involvement, a set of PARE (Parallel Analysis of RNA Ends) libraries totaling over two billion reads was constructed and sequenced from Brachypodium, switchgrass, and sorghum representing the first report of RNA degradome data from the latter two species. Analysis of this data provided not only PARE evidence for miRNA guided cleavage of over 7000 predicted target mRNAs in Brachypodium, but also evidence for miRNA guided cleavage of over 1000 homologous transcripts in sorghum and switchgrass. A pipeline was constructed to compare RNA-seq and PARE data made from Brachypodium plants exposed to various abiotic stress conditions. This resulted in the identification of 44 miRNA targets which exhibit stress regulated cleavage. Time course experiments were performed to reveal the relationship between miR393ab, miR169a, miR394ab, and their respective targets throughout the first 36 h of the cold stress response in Brachypodium. CONCLUSIONS: Knowledge gained from this study provides considerable insight into the RNA degradomes and the breadth of miRNA target conservation among these three species. Additionally, associations of a number of miRNAs and target mRNAs with the stress responses have been revealed which could aid in the development of stress tolerant transgenic crops.

An miRNA Expression Signature for the Human Colonic Stem Cell Niche Distinguishes Malignant from Normal Epithelia.

Malignant transformation of tissue stem cells (SC) may be the root of most cancer. Accordingly, we identified miRNA expression patterns in the normal human colonic SC niche to understand how cancer stem cells (CSC) may arise. In profiling miRNA expression in SC-enriched crypt subsections isolated from fresh, normal surgical specimens, we identified 16 miRNAs that were differentially expressed in the crypt bottom, creating an SC signature for normal colonic epithelia (NCE). A parallel analysis of colorectal cancer tissues showed differential expression of 83 miRNAs relative to NCE. Within the 16 miRNA signature for the normal SC niche, we found that miR-206, miR-007-3, and miR-23b individually could distinguish colorectal cancer from NCE. Notably, miR-23b, which was increased in colorectal cancer, was predicted to target the SC-expressed G protein-coupled receptor LGR5. Cell biology investigations showed that miR-23b regulated CSC phenotypes globally at the level of proliferation, cell cycle, self-renewal, epithelial-mesenchymal transition, invasion, and resistance to the colorectal cancer chemotherapeutic agent 5-fluorouracil. In mechanistic experiments, we found that miR-23b decreased LGR5 expression and increased ALDH+ CSCs. CSC analyses confirmed that levels of LGR5 and miR-23b are inversely correlated in ALDH+ CSCs and that distinct subpopulations of LGR5+ and ALDH+ CSCs exist. Overall, our results define a critical function for miR-23b, which, by targeting LGR5, contributes to overpopulation of ALDH+ CSCs and colorectal cancer. Cancer Res; 77(14); 3778-90. ©2017 AACR.

Sample sequencing of vascular plants demonstrates widespread conservation and divergence of microRNAs.

Small RNAs are pivotal regulators of gene expression that guide transcriptional and post-transcriptional silencing mechanisms in eukaryotes, including plants. Here we report a comprehensive atlas of sRNA and miRNA from 3 species of algae and 31 representative species across vascular plants, including non-model plants. We sequence and quantify sRNAs from 99 different tissues or treatments across species, resulting in a data set of over 132 million distinct sequences. Using miRBase mature sequences as a reference, we identify the miRNA sequences present in these libraries. We apply diverse profiling methods to examine critical sRNA and miRNA features, such as size distribution, tissue-specific regulation and sequence conservation between species, as well as to predict putative new miRNA sequences. We also develop database resources, computational analysis tools and a dedicated website, http://smallrna.udel.edu/. This study provides new insights on plant sRNAs and miRNAs, and a foundation for future studies.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon.

BACKGROUND: The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon. RESULTS: B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner. CONCLUSIONS: B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

MAJOR DEVELOPMENTAL REGULATORS AND THEIR EXPRESSION IN TWO CLOSELY RELATED SPECIES OF PORPHYRA (RHODOPHYTA)(1).

Little is known about the genetic and biochemical mechanisms that underlie red algal development, for example, why the group failed to evolve complex parenchyma and tissue differentiation. Here we examined expressed sequence tag (EST) data from two closely related species, Porphyra umbilicalis (L.) J. Agardh and P. purpurea (Roth) C. Agardh, for conserved developmental regulators known from model eukaryotes, and their expression levels in several developmental stages. Genes for most major developmental families were present, including MADS-box and homeodomain (HD) proteins, SNF2 chromatin-remodelers, and proteins involved in sRNA biogenesis. Some of these genes displayed altered expression correlating with different life history stages or cell types. Notably, two ESTs encoding HD proteins showed eightfold higher expression in the P. purpurea sporophyte (conchocelis) than in the gametophyte (blade), whereas two MADS domain-containing paralogs showed significantly different patterns of expression in the conchocelis and blade respectively. These developmental gene families do not appear to have undergone the kinds of dramatic expansions in copy number found in multicellular land plants and animals, which are important for regulating developmental processes in those groups. Analyses of small RNAs did not validate the presence of miRNAs, but homologs of Argonaute were present. In general, it appears that red algae began with a similar molecular toolkit for directing development as did other multicellular eukaryotes, but probably evolved altered roles for many key proteins, as well as novel mechanisms yet to be discovered.

Methods for isolation of total RNA to recover miRNAs and other small RNAs from diverse species.

For the experimental analysis of miRNAs and other small RNAs in the 20-25 nucleotide (nt) size range, the first and most important step is the isolation of high-quality total RNA. Because RNA degradation products can mask or dilute the presence of true miRNAs, it is important when choosing a method that it efficiently extracts RNA from tissues in a manner that prevents degradation of RNA of both high and low molecular weight. In addition, the presence of polyphenols, polysaccharides, and secondary metabolites may render nucleic acids insoluble, and hinder the recovery of the miRNAs. Finally, and most importantly, the method chosen must be capable of retaining the small RNA component. In this chapter, we will present a set of total RNA isolation methods that can be used to maximize the recovery of high-quality RNA to be used in miRNA analysis for a large number of plant species and tissue types.

Distinct extremely abundant siRNAs associated with cosuppression in petunia.

Cosuppression is a classical form of eukaryotic post-transcriptional gene silencing. It was first reported in transgenic petunia, where a sense transgene meant to overexpress the host Chalcone Synthase-A (CHS-A) gene caused the degradation of the homologous transcripts and the loss of flower pigmentation. In this work, we used deep sequencing technology to characterize in detail the small RNA population generated from the CHS-A sequence in cosuppressed transgenic petunia. Unexpectedly, two distinct small interfering RNAs (siRNAs) were found to vastly predominate. Our demonstration that they guide prominent cleavage events in CHS-A mRNA provides compelling and previously lacking evidence of a causative association between induction of individual siRNAs and an example of cosuppression. The preferential accumulation of these siRNAs provides new insights about sense cosuppression that may apply to other natural and engineered RNA silencing events.

Distinct size distribution of endogeneous siRNAs in maize: Evidence from deep sequencing in the mop1-1 mutant.

Small RNAs from plants are known to be highly complex and abundant, with this complexity proportional to genome size. Most endogenous siRNAs in Arabidopsis are dependent on RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) for their biogenesis. Recent work has demonstrated that the maize MEDIATOR OF PARAMUTATION1 (mop1) gene is a predicted ortholog of RDR2. The mop1 gene is required for establishment of paramutation and maintenance of transcriptional silencing of transposons and transgenes, suggesting the potential involvement of small RNAs. We analyzed small RNAs in wild-type maize and in the isogenic mop1-1 loss-of-function mutant by using Illumina's sequencing-by-synthesis (SBS) technology, which allowed us to characterize the complement of maize small RNAs to considerable depth. Similar to rdr2 in Arabidopsis, in mop1-1, the 24-nucleotide (nt) endogenous heterochromatic short-interfering siRNAs were dramatically reduced, resulting in an enrichment of miRNAs and transacting siRNAs. In contrast to the Arabidopsis rdr2 mutant, the mop1-1 plants retained a highly abundant heterochromatic approximately 22-nt class of small RNAs, suggesting a second mechanism for heterochromatic siRNA production. The enrichment of miRNAs and loss of 24-nt heterochromatic siRNAs in mop1-1 should be advantageous for miRNA discovery as the maize genome becomes more fully sequenced.